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O ' ABSTRACT 

, We address the question of whether the angular power spectrum of the thermal 

' Sunyaev-Zel'dovich (SZ) sky is further distorted by weak gravitational lensing of fore- 

O , ground large-scale structures. Using an analytic approach to both gaseous and dark 

' halo models, we show that the contamination of weak lensing in the measurement of 

SZ power power is negligibly small, and relatively corrections |AC;|/C; are less than 
OO ' 3% up to ^ = 10^. This arises from both the weaker gravitational potentials of low- 

redshift matter inhomogeneities that can act as lenses for SZ sources (clusters) and 
the shallower shape of intrinsic SZ power spectrum at large I, in contrast to the cosmic 
microwave background which can be significantly affected by weak lensing because of 
Q>^ ' the distant location and significant damping of its intrinsic power spectrum at small 

\Q ', angular scales. 

■ Key words: cosmic microwave background — galaxies: clusters: general - gravita- 

tional lensing — intergalactic medium — large-scale structure of universe 
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Q_i 1 INTRODUCTION 

O Gravitational lensing by foreground large-scale structures of the universe distorts many varieties of background sources. Among 

I these, detection of weak lensing-induced signatures in the cosmic microwave background (CMB) will provide a sensitive probe 

52 ' of matter distribution on large-scales and at high redshifts up to z « 1100 (e.g. Zaldarriaga & Seljak 1999). On the other 



> 



hand, the CMB photons can also be distorted by non-gravitational effect along the path from the last scattering surface 
to us. The most significant pattern is the excess of the CMB power spectrum relative to primordial CMB anisotropy at 
I > 2000 (Dawson et al. 2001; Bond et al. 2003), which is produced by the inverse Compton scattering of the CBM photons 



?H ' by intervening energetic electrons confined within massive clusters, known as the Sunyaev-Zel'dovich (SZ) effect (Sunyaev & 
^ Zeldovich 1972). Studies of both the SZ-selected clusters and the SZ power spectrum yield important information about the 
distribution and evolution of hot intracluster medium in the universe. 

At scales below 10 arcminutes characterized hy I ~ 2000-10*, the CMB power spectrum is dominated by the thermal 
SZ effect of clusters distributed at intermediate/high redshifts up to 2 « 1 (Cooray 2000; Cooray & Ifu 2000). An interesting 
question thus arises: Can the SZ sky be further distorted by the weak gravitational lensing effect of foreground matter 
inhomogeneities ? Namely, we would like to address the question of to what extent the SZ power spectrum at small angular 
scales should be corrected for because of weak lensing by large-scale structures. It has been noted that the deep gravitational 
potentials of clusters can lead to both steplike wiggles on the CMB sky (Seljak & Zaldarriaga 2000) and a systematic 
magnification bias of unresolved background radios behind clusters (Loeb & Refregier 1997). The latter may contaminate the 
measurement of the thermal SZ effect if the removal of background radio sources is made without inclusion of the inevitable 
lensing magnification bias. Furthermore, the SZ sky is also correlated with the weak lensing map because the two effects are 
generated by common massive dark halos (Goldberg & Spergel 1999; Cooray & Hu 2000). The problem we explore in this 
paper is more straightforward, in the sense that the background source is now the SZ map which is associated with each cluster, 
while the foreground lenses are large-scale structures along the line of sight. Technically, we can use the formalism similar to 
the calculation of weak lensing of the CMB power spectrum, which has been well developed in the literature (Seljak 1996; 
Bernardeau 1997; Zaldarriaga & Seljak 1999; Hu 2000; Zaldarriaga 2000), except that the distribution function of background 
sources (clusters) should now be taken into account for our case. Throughout this paper we assume a fiat cosmological model 
(ACDM) with the parameters determined recently by WMAP: Q,m = 0.27, SIa = 0.73, flth'^ = 0.0224, h = 0.71, Us = 0.93 
and as = 0.84. 
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2 FORMALISM 



Following the treatment of weak lensing of the CMB under the flat-sky approximation (e.g. Bernardeau 1997; Zaldarriaga & 
Seljak 1999; Hu 2000; Zaldarriaga 2000), we calculate the lensed SZ temperature fluctuation T{0) in the direction 6 by 



oc 

T{e) = r''(0 + v</-) = ^ i [V0(6») • vYT%e), 



(1) 



where T"(0) is the unlensed SZ temperature pattern in direction produced by clusters at redshift Zs, V(j> denotes the lensing 
deflection angle by foreground large-scale structures along the line of sight, and the projected potential is related to the 
three-dimensional potential $(r) of density perturbation S{r) at z through 
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in which Dds and Ds are the angular diameter distances from the matter inhomogeneities at z and the observer at z = to the 

sources at Zs, respectively, r is the comoving distance and = (1 -I- z). The power spectrum of the lensed SZ temperature 
map up to the fourth order in deflection angle in terms of equation (1) is thus (e.g. Cooray 2003) 



Ci 



Cf - Gi 



J i^^y 



■II 



2 |1'|4 



I 



(277)2 

dH' 



dH" 



■IT 



ii'i* 



(27r)2 (27r)2"^l'-''l^''^'" 
[(1 - 1' - 1") ■ If [(1 - 1' ■ 



'^I'^l" 



(1 - 1') ■ If [(1 - 1') • I'f 



1" 



|1'|4 



|1„|4 



2/ (27i-)2 (27r)2 I*-* ' I |1'|4 |1"|4 

where Cf is the unlensed power spectrum of the SZ map, and Cf corresponds to the lensing power spectrum defined as 



Cf = 



s)dZs 
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(4) 
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10 Jo 

«dec ~ 1100 is the decoupling redshift. Da is the comoving angular diameter distance, P™(fc, «) is the matter power spectrum 
at z, and 'w{zs) is the distribution function of the SZ sources (i.e. clusters). If the mass function of dark halos is d^n/dMdV, 
then w{zs) can be expressed as 



dzsdfi J M. 
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(6) 



where d^V / dZgdO, denotes the comoving volume per unit redshift and per steradian. We have also tested the observationally 

determined X-ray luminosity function of clusters (Rosati, Borgani & Norman 2002) for Lx > 10*'' erg s~^ instead of the mass 
function d^n/dMdV and found that the modification is only minor. In principle, the lower mass limit Mmin can be taken to 
be because low- mass halos (e.g. M < IQ^^Mq) almost make no contributions to SZ effect unless the self-similar model is 
adopted for gas distribution. Nevertheless, inclusion of the low-mass halos may significantly alter the distribution function 
w{zs) ■ So, we will take Mmin = 10^** Mq in the following calculation, which corresponds roughly to a lower luminosity cutoff 
of Lx > W*^ erg s~^. Such a choice of the lower limit is also consistent with the prediction of preheating model, in which the 
preheated baryons cannot be trapped in the shallower gravitational potential wells of low-mass halos (sec below). 

The undistorted thermal SZ power spectrum Cf is composed of the so-called Poisson term Cf and clustering term Cf , 
respectively. 
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(8) 



where yi{M,Zs) is the Fourier transform of the Compton j/-parameter of the cluster of mass M at Zs, -Do = (1 + Zs)Ds 
is the comoving distance, P''°(fc, 2s) is the linear matter power spectrum at Zs, and h{M,Zs) is the bias parameter, for 
which we use the analytic prescription of (Mo & White 1996). The frequency dependence term g{x) is defined as g{x) — 
(a;^e^)[4 — a;coth(a;/2)]/(e^ — 1)^, in which x = hu/kTcMB is the dimensionless frequency, and Tcmb is the temperature of 
the present CMB. The observing frequency will be fixed at u = 30 GHz in this study. 

Mass function of massive dark halos is assumed to follow the standard Press-Schechter function 
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W{kR) is the Fourier representation of the window function, and W{kR) = 3([sin(feJ?) — kRcos{kR)]/ (kR)^ . We adopt the 
present-day linear matter power spectrum P''"(fe) given by Baxdeen et al. (1986) and a power-law form (oc fe"") for the 
primordial matter density fluctuation. 

At small scales below ^ 10', matter power spectrum is strongly affected by non-linear structures. Following Scljak (2000) 
and Cooray, Hu & Miralda-Escude (2000) we use a halo approach to calculate the non-linear dark matter power spectrum 
P™(fe, z), which consists of a single halo term P^^{k, z) and a clustering term P'^^{k, z): 



where Ph{M, z, k) is the Fourier transform of dark matter density profile ph{M, z) of a halo M at z, and p{z) is the mean mass 
density of the background universe at z. We adopt the functional form suggested by numerical simulations (Navarro, Frenk 
& White 1997; NFW) for ph{M,z) and fix the concentration parameter in terms of the empirical fitting formula of Bullock 
et al. (2001). 

Finally, we need to specify the gas distribution of dark halos in order to evaluate the Compton y-parameter. Since the 
weak lensing effect of the SZ map may occur at very small angular scales, wc should employ a more realistic model for the gas 
distribution in the central regions of clusters instead of the commonly adopted self-similar model in literature. Recall that the 
SZ power spectra predicted by the self-similar model and other physically-motivated models such as preheating and radiative 
cooling become indistinguishable roughly at I < 2000. For simplicity, here wc work with an analytic preheating model, in 
which the gas distribution in massive halos can be obtained by combining the equation of hydrostatic equilibrium with a 
modified entropy profile S{r) = Sfloor + iSNFw(r'). Namely, a constant entropy floor of Sfloor ^ 120 keV cm^ is now added 
to the entropy profile S'nfw('^') predicted by the self-similar model which again takes the universal density profile (NFW) . 
This model allows us to essentially reproduce the observed X-ray surface brightness distribution (i.e. (3 model) and statistical 
properties of groups and clusters (e.g. the X-ray Luminosity - temperature relation and entropy distribution) (Babul et al. 
2002; Xue & Wu 2003). 



3 RESULTS 

We first demonstrate in Figure 1 the lensing power spectrum Cf for the SZ sky. For comparison, we have also shown the 

result for the CMB calculated by setting w{zs) — 5{zs — 2;doc) in eqn.(5) and running the publicly available code CMBFAST. 
At small scales, the lensing power spectra for the SZ and CMB maps appear to be similar in shape, but the lensing amplitude 
for the SZ sources is about two orders of magnitude smaller than that for the CMB. This is simply because there are much 
more intervening matter inhomogeneities acting as lenses for the CMB photons than those for the SZ signals generated by 
clusters located mostly at 2; < 1. Furthermore, Cf peaks at lower I for the SZ map than for the CMB, consistent with the fact 
that the large-scale structures responsible for weak lensing of the SZ sky are at low redshifts and therefore, subtend relatively 
large angles. Recall that gravitational lensing conserves the surface brightness of background source. Therefore, Cf acts as 
a narrow window function on the background power spectrum Cf via eqn.(4), which alters the intrinsic patterns of Cf . Our 
intuition in terms of this simple argument is that the weak lensing effect on the SZ power spectrum is probably much smaller 
than the lensing signature on the CMB, if the lensing induced correction is simply proportional to Cf . Finally, the lensing 
power spectra at small scales are dominated by non-linear structures, and linear approximation of matter power spectrum 
P''°(A;, z) may gives rise to an underestimate of the lensing effect. 

Nevertheless, the extent to which background sources are affected by weak lensing of foreground large-scale structures 
also depends critically on the shape of the source power spectrum. We begin by recalling how the CMB can effectively be 
lensed by large-scale structures (see Hu 2000): At large scales I < 100, the CMB power spectrum can be considered to be 
a slowing varying function of I (see Figure. 2), and we can approximately take Cf out of the second integral in eqn(4). The 
lensing induced two terms to the first order approximation in defiection angle cancel if 1 — 1' is further replaced by 1. This 
yields Ci « Cf . At small scales, existence of the strong Silk damping of intrinsic CMB power spectrum at I > 2000 makes the 
contribution of the first correction term in eqn.(4) considerably small, and the second term dominates the lensing correction. 
As a result, weak lensing always leads to an increase of the CMB power at small angular scales, as shown in Figure 2. Now 
we turn to the SZ power spectrum: Unlike the intrinsic CMB, there is no significant damping effect on the thermal SZ power 
spectrum over all angular scales. Rather, the SZ power spectrum has a long, slowly decreasing tail at high I. In other word, the 
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Figure 1. Power spectra of weak lensing for CMB and SZ map. Dotted and solid lines represent the results of linear and non- linear 
matter power spectra, respectively. 



Table 1. Numerical Results 



I 


100 


1000 


10000 


100000 


"^'"5' (%) 

I 

(fourth order) 
Ci (second order) 


0.0025 
7.9 X 10-10 


0.0065 
9.2 X 10-» 


0.53 
1.5 X 10-5 


2.51 
6.7 X 10-" 



first negative term in eqn.(4) can no longer be ignored. Applying the approximation that Cf is slowly varying and 1 — 1' « 1 
to eqn.(4), we have C; « C" to the first order in deflection angle. Our numerical calculation reaches essentially the same 
conclusion. In Table 1 we have listed the relative corrections \Ci — Ci\/Ci for a sot of multipoles from I = 100 to I = 10^. It 
turns out that the lensing induced corrections are well below 3% even up to Z « 10^. Meanwhile, we have given in Table 1 
the power spectrum ratios of the fourth order correction [the last three terms of eqn.(4)] to the second order approximation 
in cf [the second and third terms of eqn.(4)[. While the eoutribution of higher order corrections to the power spectrum Ci 
becomes increasingly important with the decrease of angular scale, the effect of higher order terms on the evaluation of Cj is 
still negligibly small on scales up to / = 10® as compared with the contribution from the first order approximation in defiection 
angle. 



4 CONCLUSIONS 

It has been shown that the angular power spectrum of the thermal SZ sky is almost unaffected by the weak lensing of 

foreground large-scale structures. This is primarily due to both the relatively weaker gravitational potentials of low-redshift 
matter inhomogeneities that can act as lenses for SZ sources and the shallower shape of SZ power spectrum at smaller angular 
scales. Note that we have adopted the gas density profile predicted by the so-called preheating model, which results in a more 
rapid decline of the SZ power spectrum at large / than the one predicted by the self-similar model or cooling model (Wu & Xue 
2003). In other word, we have already provided an optimistic estimate of the weak lensing effect on the SZ power spectrum 
since significant fluctuations intrinsic to background sources are more easily be distorted by lensing. This explains why weak 
lensing of large-scale structures of the universe can effectively alter the CMB sky at Zs 1100 as the CMB power spectrum 
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Figure 2. Angular power spectra of the unlensed and lensed CMB and SZ map. The lensed SZ power spectrum is almost identical to 
the intrinsic one, and the relative difference rea<;hes only 3% at I fs 10^. 



exhibits the sharp damping signatures at high I (Hu 2000). The angular power spectrum of 21cm map of the reionization at 
« ~ 10 also demonstrates a significant decUne feature at I = 10^ — 10'' (e.g. Zaldarriaga, Furlanetto & Hernquist 2003), and 
it is thus expected that the 21cm fluctuation of the reionization should be affected by weak lensing (Pen 2003). Nonetheless, 
the present study has eliminated the concern that the SZ sky may be further contaminated by weaJj lensing of foreground 
large-scale structures of the universe. 
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